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The effect of a weak uniform magnetic field on the electronic structure of slightly deformed 
fullerene molecules is studied within the continuum field-theory model. It is shown how the exist- 
ing due to spheroidal deformation fine structure of the electronic energy spectrum modifies in the 
presence of the magnetic field. Exact analytical solutions for zero-energy modes are found. 
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I. INTRODUCTION 



Recently, we have considered the problem of the low energy electronic states in spheroidal fullerenes [1]. The main 
findings were a discovery of fine structure with a specific shift of the electronic levels upwards due to spheroidal 
deformation. In addition, three twofold degenerate modes near the Fermi level with one of them being the true 
zero mode were found. An interesting question is the modification of this structure under the influence of a uniform 
magnetic field. 

| The problem of Zeeman splitting and Landau quantization of electrons on a sphere was studied in Ref. [2]. To 
c"j . this end, the Schrodinger equation for a free electron on the surface of a sphere in a uniform magnetic field was 
formulated and solved. In this paper, we explore the field-theory model suggested in Ref. [1], which describes the 
electronic states near the Fermi energy and takes into account the specific structure of carbon lattice, geometry, and 
the topological defects (pentagons). The Euler's theorem for graphene requires the presence of twelve pentagons to 
get the closed molecule. In the framework of continuum description we extend the Dirac operator by introducing the 
t> ■ Dirac monopole field inside the spheroid to simulate the elastic vortices due to twelve pentagonal defects. The K spin 
flux which describes the exchange of two different Dirac spinors in the presence of a conical singularity is included in 
a form of t'Hooft-Polyakov monopole. Our studies cover slightly elliptically deformed molecules in the weak uniform 
external magnetic field. 
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O , II. THE MODEL 



Let us start with writing down the Dirac operator for free massless fermions on the Riemannian spheroid S 2 . The 
Dirac equation on a surface £ in the presence of the abelian magnetic monopole field and the external magnetic 
field is written as [3] 

i'fe^\S7 lt -iW fl -iA lt ]i> = Ei>, (I) 

where is the zweibein, = e^e^Sap is the metric, the orthonormal frame indices a, (3 = {1,2}, the coordinate 
indices /i, v = {1,2}, and V M = + with 

% = 1^^,713}, (2) 

being the spin connection term in the spinor representation (see [I] for details). The energy in (I) is measured from 
the Fermi level. 



A. Zero-energy mode 

Let us start from the analysis of an electron state at the Fermi level (so-called zero-energy mode) . We will use the 
projection coordinates in the form 
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2R 2 r , 2R 2 r , R 2 - r 2 

where R and c are the spheroidal axles. The Riemannian connection with respect to the orthonormal frame is written 

as [4,5] 



R 2 -r 



2 



^2 - ">1 - „ 2 v 2 , .on ' "V2 - "rl 



0. (4) 



v /( i? 2_ r 2)2 +4c 2 r 2 

The only nonzero component of the gauge field in region R N for spheroidal fullerenes reads (see Ref. [1]) 



where 



W^G + m^, (5) 

T^Tr 2 ' (6) 

Notice that the monopole field in Eq. (5) consists of two parts. The first one comes from the K-spin connection 
term and implies the charge g — ±3/2 while the second one is due to elastic flow through a surface. This contribution 
is topological in its origin and characterized by charge G. For total elastic flux from twelve pentagonal defects one 
has G = 1. The parameter m — g — G is introduced in Eq.(5). The uniform external magnetic field B is chosen to be 
pointed in the z direction so that A = B (y, — x,0) /2. In projection coordinates the only nonzero component of 
is written as 

2BR 4 r 2 

Let us study Eq.(l) for the electronic states at the Fermi energy (E = 0). The Dirac matrices can be chosen to be 
the Pauli matrices, 71 = — 02,72 = —o'l- By using the substitution 

^)_r e^ f^)) ,--0±l±2 (8) 

1/ z . x R 2 + r a /n r 2 -R 2 A , . n 

1 z R 2 + r 2 r 2 - R 2 \ 

-<, _ ,„_ _ a 4 ) - _+_(* - sp^)) »W = 0. 

where K = \J (R? — r 2 ) 2 + Ac 2 r 2 . We assume that the eccentricity of the spheroid is small enough. In this case, 
one can write down c = R + SR and 5 (| 6 |<C 1) is a small dimensionless parameter characterizing the spheroidal 
deformation. Therefore, one can follow the perturbation scheme using S as the perturbation parameter. To the leading 
in 6 approximation Eqs. (11) and (12) are written as 

(R 2 + r 2 ) 2 r 2 - R 2 (r 2 - R 2 ) 3 

^Ar) = (- r[{R2+r2) 2_ 4R 2 r 2 S] ij - + *)^R2 + m \ r 2 + R 2 r ] 

2BR 4 r r 2 - R 2 



we obtain 



[{R 2 + r 2 ) 2 - 4R 2 r 2 S] 2r(R? + r 2 ) 



>j(r), (11) 



* ,s , (R 2 + r 2 ) 2 r . ,., „r 2 -R 2 Jr 2 -R 2 )\ 

d r uAr) = ( 5^ [j - mil + b)^ — T + md—^ — r-^ 

Jy ' r[(R 2 + r 2 ) 2 - AR 2 r 2 S] r 2 + R 2 (r 2 + R 2 ) 3 ' 

2BR 4 r r 2 -R 2 

+ [ (i? 2 + r 2)2 _ AR 2 r 2$} + 2r (i?2 + r 2) ) U i( r J- ( 12 ) 
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The exact solution to Eqs. (11) and (12) is found to be 



^ (.x + l^/a-m) x + i _ 25 - 2^8(8 - 1) \ m 
f 7 -(x) = w . / v , ((a; + 1) - 4fe) , (13) 



where a; = r 2 /i? 2 , a = (j/2)y/5/(5 - 1) + BR 2 /^8(8 - 1). The function can be normalized if the condition 
m — 1/2 < j < 1/2 — to is fulfilled. In turn, the normalization condition for vj has the form —l/2 — m<j<m+l/2. 
Finally, there are five normalized solutions for Uj with j = 0, ±1,±2 when to = —5/2, and one solution for Vj with 
j = when to = 1/2. In the limit 8 — > we arrive at the solution for spherical fullcrcnes 

Wj -(r) = A- 1 ' 2 -?-^ (R 2 + r 2 )( 1 /2+H e ^ ; 
Uj (r) - r (-V2+j+m) (#2 + ^(i/2-m) e -^n ( 15 ) 

Since normalization conditions do not depend on the parameter 6 they are the same as for spherical case. Thus the 
perturbation does not change the number of zero modes. Notice that Eq. (15) can be obtained from Eqs. (13) and 
(14) by using lim [1 + (x/a)] a — e x , a — > oo. 



B. Electronic states near the Fermi energy 



Let us study now the structure of electron levels near the Fermi energy. It is convenient to consider this problem 
by using of the Cartesian coordinates x, y, z in the form 

x = a sin 9 cos (f>; y = a sin 9 sin <fr; z — c cos 9. (16) 

The Riemannian connection reads (cf. (4)) 



a cos 9 

= -"01 = , ?n =f= ; "92 = "81 = 0. (17) 

V a cos z 9 + cr sin 1 



Within the framework of the perturbation scheme the spin connection coefficients are written as 

^ 2 = -^i«cos0(l-(5sin 2 #), (18) 

where the terms to first order in 8 are taken into account. In spheroidal coordinates, the only nonzero component of 
Wfj, in region Rn is found to be (see Ref. [1]) 

w g cos 9(1 + 8 sin 2 9) + G(l - cos 9) - 8G sin 2 9 cos 9, (19) 

and the external magnetic field reads 

A^^-^Ba 2 sin 2 9. (20) 
By using the substitution (8) we obtain the Dirac equation for functions Uj and Vj in the form 

+ + i*> ^ ,) + *,) ( ) - E ( ) . (21, 



where 



A = - — sin 9 (j - 2m cos 9) - 71 sin 3 9. (22) 

CI A 
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A convenient way to study the eigenvalue problem is to square Eq. (21). For this purpose, let us write the Dirac 
operator in Eq. (21) as V — T>o + 8V\. One can easily obtain that (similar way as in Ref. [6]) 

^9 1/^9 cos9 n 1 1 \ (j — mcos9) 2 „ B , „ , 

a 2 \ ° sin 6 4 4 sin 2 / a 2 sin' 1 9 2 

m - j cos B 2 a 2 sin 2 . . 

+(T 3 2 7 2fl + 7 • (23) 

<r sin 9 4 

Notice that to first order in S the square of the Dirac operator is written as V 2 = {V\ + <5r), where T = (V T>i +V{D ). 
In an explicit form 

a 2 f = 2j 2 +jx (cr 3 - 6m) + 4m(m - a 3 )x 2 + 2ma 3 + 3Ba 2 (l - x 2 )x(a 3 /2 - m) 

+2Bja 2 (l - x 2 ) + B 2 a 4 {l - x 2 ) 2 /2 (24) 
where the appropriate substitution x = cos 9 is used. The equation V 2 ip = E 2 ip takes the form 

- - 2 )^ - (j "^ )2 7;:r i+g3m - « 2 ^(.) - ^ 2 r] ( ^ 

- -(A 2 - 7) f ) , (25) 



4 y V 

where A = aE, V(x) = j + (<r 3 — 2m)x/2. Since we consider the case of a weak magnetic field the terms with B 2 and 
5B in Eq. (25) can be neglected. As a result, the energy spectrum for spheroidal fullcrcncs is found to be 

(X s jn ) 2 = (n + \j\ + 1/2) 2 - to 2 + Ba 2 j + Ba 2 A,„ + <5F(j, n, m) (26) 

where 

j(m 2 - 1/4) 



A — — - 



P(P+1) 



The function of spheroidal deformations F(j, n, m) is specified in Ref. [1]. Finally, in the linear in 8 approximation, 
the low energy electronic spectrum of spheroidal fullercnes takes the form 



with 



E jn - ±V (24 + n)(2r] + n), E jn - , E jn , 

jn jn 

where 4 = A* i v ) an( i V — ft ( a ) f° r 3 > (j < 0), respectively. Here we came back to the energy variable E = X/a (in 
units of KVf/cl where Vp is the Fermi velocity). 

As is seen from Eq. (27), for B = the spheroidal deformation gives rise to an appearance of a additional 
structure. As an example, Table 1 shows all five contributions (in compliance with Eq. (27)) to the first and second 
energy levels for YO-C240 (YO means a structure given in [7,8]). As is seen, the energy levels become shifted due 
to spheroidal deformation. The uniform magnetic field provides the well-known Zeeman splitting. The difference 
between topological and Zeeman splitting is clearly seen. In the second case, the splitted levels are shifted in opposite 
directions while for topological splitting the shift is always positive. As an illustration, Fig. 1 schematically shows 
the structure of the second level. In this case, the initial (for (5 = 0) degeneracy of E® n is equal to six. The spheroidal 
deformation provokes an appearance of three shifted double degenerate levels (fine structure). The magnetic field is 
responsible to Zeeman splitting. 

We admit that the obtained values of the splitting energies can deviate from the estimations within some more 
precise microscopic lattice models and density- functional methods (see, e.g. [9]). In our paper, we focused mostly on 
the very existence of physically interesting effects. For this reason, the values in both Table and the schematic picture 
are presented with accuracy at about one percent of Ej n . 
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III. CONCLUSION 



In conclusion, we have considered the structure of low energy electronic states of spheroidal fullerenes in the weak 
uniform magnetic field provided the spheroidal deformation from the sphere is small enough. For the states at the 
Fermi level, we found an exact solution for the wave functions. It is shown that the external magnetic field modifies 
the density of electronic states and does not change the number of zero modes. For non-zero energy modes, electronic 
states near the Fermi energy of spheroidal fullcrcne arc found to be splitted in the presence of a weak uniform magnetic 
field. It should be mentioned that the zero-energy states in our model corresponds to the HOMO (highest occupied 
molecular orbital) states in the calculations based on the local-density approximation in the density-functional theory 
(see, e.g. [10]). In particular, the HOMO-LUMO energy gap is found to be about 1.1 eV for YO-C240 fullerene within 
our model (here LUMO means the lowest unoccupied molecular orbital). 

The work was supported in part by VEGA grant 2/7056/27. of the Slovak Academy of Sciences, by the Science 
and Technology Assistance Agency under contract No. APVT-51-027904 and by the Russian Foundation for Basic 
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YO - C240 Ba 2 = 0.1 


J 


E%{eV) 


E™°(meV) 


E$ n (meV) 


n = 0, m = 1/2 


1 

-1 


1.094 
1.094 


27 
-27 


10.5 
10.5 


n = 1, m = 1/2 
n = 0, m = 1/2 
n = 0, m = —5/2 


1 

-1 
2 
-2 
3 
-3 


1.89 
1.89 
1.89 
1.89 
1.89 
1.89 


16 
-16 
32 
-32 
24 
-24 


8.8 
8.8 
28.4 
28.4 

3 

3 



TABLE I. The structure of the first and higher energy levels for YO-C240 fullerene in uniform magnetic field. The hopping integral 
and other parameters are taken to be t = 2.5 eV and V F = 3tb/2h, b = 1.45A, R = 7.03A, SD = 0.17A, 8 = 0.024. b is average bond 
length, R (R = a) is average radius, SD is standard deviation from a perfect sphere (see Refs. [7,8]), so that 5 = SD/R. 



E% +E$ n +E° 
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FIG. 1. The schematic picture of the second positive electronic level Ej n of the spheroidal fullerenes in a weak uniform 
magnetic field. 
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